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Abstract 

We study unparticle effects on particle and antiparticle osillations in meson-antimeson, and 
muonium-antimuonium systems. Unlike usual tree level contributions to meson oscillations from 
heavy particle exchange with small Ti2, the unparticle may have sizeable contributions to both Myi 
and Ti2 due to fractional dimension du of the unparticle. We find that very stringent constraints 
on the unparticle and particle interactions can be obtained. If unparticle effect dominates the 
contributions (which may happen in D° — D° mixing) to meson mixing parameters x and y, we 
find that x/y = cot(irdu). Interesting constraints on unparticle and particle interactions can also 
be obtained using muonion and antimuonion oscillation data. We also comment on unparticle 
effects on CP violation in meson oscillations. 



I. INTRODUCTION 



Recently Georgi proposed an interesting idea to describe possible scale invariant effect at 
low energies by unparticlesjX]. Georgi argued that operators Obz made of BZ fields in the 
scale invariant sector may interact with operators 0$m of dimension d$M made of Standard 
Model (SM) fields at some high energy scale by exchange particles with large masses, My, 
with the generic form OsmObz/M^. At another scale the BZ sector induce dimensional 
transmutation, below that scale the BZ operator Obz matches onto unparticle operator Ou 
with dimension du and the unparticle interaction with SM particles at low energy has the 
form 

X^- dsM - du Os M O u . (1) 



Study of unparticle effects has drawn a lot of attentions 



rom collider physics 



% |5|], long range effects 



2], 



6], 



low energy flavor conserving and flavor violating processes j^l 

cosmological and astrophysics phenomena (7|, to more theoretical*,! hi. .i™* 

we further study unparticle effects on particle and antiparticle oscillations of meson P and 
antimeson P, and muonium and antimuonium systems. 

Unparticle effects on oscillation in meson and antimeson have been considered previ- 
ously [4]. Our investigation for meson and antimeson oscillation will focus on some inter- 
esting features due to the fractional dimension of unparticle du- Unlike usual tree level 
contributions to meson oscillations from heavy particle exchange with small T 12 , the un- 
particle may have sizeable contributions to both M\2 and due to fractional dimension 
du of the unparticle leading to a phase factor (— \) du ~ 2 in the propagation. If unparticle 
effect dominates the contributions (which may happen in D° — D° mixing) to meson mixing 
parameters x and y, we find that x/y = cot(irdu)- 

Meson- antimeson oscillation exists in several neutral meson systems, K° — K°, D° — D°, 
B® — B® and B° s — B° s . Long distance contributions to oscillation parameters for K° — K° 
are large which causes a large uncertainty in theoretical calculations. We will restrict our 
calculations for unparticle effects to D° — D° and B® s — B® s systems. We find that very 
stringent constraints on the unparticle and particle interactions can be obtained. 

Muonium (M = (p,e)) and antimuonium (M = (e/x)) oscillation may also provide inter- 
esting constraints on flavor changing interaction. Experimentally muonium-antimuonium 
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oscillation has not been established. Our analysis shows that constraints on unparticle 
and particle interactions can indeed be obtained using experimental data on muonium- 
antimuonium oscillation. 

II. MESON AND ANTIMESON OSCILLATIONS 

The mixing of a meson and its antimeson is determined by the off diagonal matrix ele- 
ments M\2 and Ti2 in the Hamiltonian. Their relations to the mass and lifetime differences 
are given by, 

K - m L ) - i{T H - Y L )/2 = 2y / (M 12 - iT l2 /2){M* 12 - iY\ 2 /2\ (2) 

where the subscripts "H" and "L" label the mass eigenstates, \Ph) = p\P) + Q\P) an d 
\Pl) — p\P) ~ Q.\P)i respectively, p and q are normalized as \p\ 2 + \q\ 2 = 1 and {q/p) 2 = 
(M* 2 — iT\ 2 /2) / (M12 — iT 12/ 2). We denote the mass and lifetime differences by Am = 
niu — m L and Ar = Th — T^. The parameters x and y are related to Am and Ar by 
x = Am/r and y = Ar/2I\ 

There are several possible contributions to M\2 and Ti2 from unparticle and particle 
interactions. The following operators composed of SM fields and derivatives with dimensions 
less than or equal to 4 invariant under the SM gauge can contribute to meson mixing at tree 
level, 

a) : X'qqAIT^Ql^QlO^ X' uv A^U r1i JJ r O^ \ ! DD ^ D Rl ,D R 0^ 

b) : iX QQ A u d «Q Lltl D»Q L U) iXuuAu du URl^U R Ou, iX DD A u du D Rl ,D»D R O u ; 

c) : t~X QQ A^Q Llfl Q L d^O u , fXuuAu^URl^Rd^Ou, i\ DD A u du D Rl ,D R d»O u ; 

d) : X YU A u d "Q L HU R O u , X YD A u d «Q L HD R O u . (3) 

Here Ql, U r , and D R are the SM left-handed quark doublet, right-handed up-quark, and 
right-handed down-quark, respectively. 

After using equation of motion for quarks, the interactions in eq. ([3]) can be parameterized 
in the following form 

For vector 0% : £ v = A l u dli q^c^^L + c VR l^R)qjO^ ; 

For scalar O u : C s = A^^-q^L + 4 R R) qj O^ , (4) 
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FIG. 1: The t and s channel contributions to meson-antimeson oscillation. 

where q^j) denote quarks with flavor indices i(J). The parameters c Vl and Cg are different 
for the interactions listed in a) to d). They are given by 



a) : 


C V L 


ij 

— A QQi C V R — A UU,DD) 
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d) : 




\ ji* ij \ ij 

= v\y F , c4 = v\J F . 









(5) 

where v = (H) is the vacuum expectation value of H. Note that the vector unparticle 
couplings cy scaled as u , while scalar unparticle couplings cs scaled as K7 11 
Evaluating the two diagrams in Fig. [TJ we obtain 
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H eff 
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H eff 



A, 



2 / 
« *2(l-d M ) m j -iirdu 1 ( ]_ 

ITT W TT e 7 7371 
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^7) Ws l l + 4 r RV 



qA . (6) 



2sin(7rd w )~~" A&~ 4 V(s) 2 " d " (t) 2 " 

Here A dw = (167r 5 / 2 /(27r) 2dw )r(d w + l/2)/(T(d u - l)T(2d u )). We have used 
(zA dw /2sin(7r^)) x (l/(-p 2 f- d u) and (^ dw /2 sin(vr^)) x {{-g^ + jfp" / J?) / {-p 2 ) 2 -"") 
for scalar and vector unparticle propagators, respectively. 

In the systems we are studying, mesons are made of a light (labelled by i) and a heavy 
quark (labelled by j). In the heavy quark limit, one has s = t « m 2 w m|>. With this 
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approximation and theoretical matrix elements for the relevant operators, we have 

ry 2 = 2Mgtan(c!„7r). (7) 



We have included a missing factor of 1/2! due to Wick rotation in previous studies 
The parameters By,s are the bag factors which are equal to 1 in the vacuum saturation and 
factorization approximation. 

We would like to point out some silent features of the unparticle contribution to M^ 2 
and r^ 2 due to the phase factor e~ mdu . We note that M^ 2 can have both sign depending 
on the value of du due to the factor cot(7r<i^), therefore if information about the sign can 
be obtained from other considerations, the dimension du can be restricted. There may be 
a sizeable contribution to Ti2 at tree level which is not possible for usual tree level heavy 
particle exchange. For du equal to half integers, there is no contribution to Mi 2 , but there is 
for T 12 . Another interesting feature of unparticle contribution is that the ratio M 12 /(r 12 /2) 
of unparticle contribution is related to the unparticle dimension parameter du by 

'- --cot(vrG^). (8) 



If the unparticle contribution dominates meson and antimeson oscillation then the mea- 
surements of M\2 and Ti 2 provide a possible way to determine the dimension parameter 
d u - 

We now present our numerical results on the constraints for unparticle and particle in- 
teractions with the assumption of CP conservation. In this case the unparticle contribution 
to mass difference Am u = 2|M^| and x u /y u = M^/(r^ 2 /2) = cot(ird u ). We will comment 
on possible CP violating effects later. The results are shown in Figs. [2] and [31 The con- 
straints are obtained with the vacuum saturation approximation, i.e. By = Bs = 1 and 
f D = 0.201 GeV, f Bd = 0.216 GeV and f Ba = 0.260 GeV [9|. 
D° - D° System 

Belle and BABAR collaborations have recently published evidence for D° — D° oscil- 
lation 10]. The Heavy Flavor Averaging Group (HFAG) [11] combined all mixing mea- 



surements to obtain world average (WA) values for x and y for CP conserving case with 
x = (0.87to;34)% and y = (0.66lo;2o)%. Short distance contributions in the SM are several 



orders of magnitudes smaller than the central experimental values. There are possible large 
long distance contributions which are however difficult to have precise predictions. We will 
assume that the contributions to x and y are purely from unparticle effects. With this 
assumption, we immediately obtain 

- = cot(7rd w ) ~ 1.31 ± 0.61 . (9) 

Although the sign of x and y are both positive, the absolute sign of cannot be 
determined from x measurement, therefore ndu can be in the first and third quadrants from 
the sign of x/y. We have 

d u = (0.21 + n) ±0.07, (10) 

where n is an integer number which cannot be determined from just information from x/y. 
Note that the experimental errors are still large, consequently the uncertainties of d u are 
substantial. 

One can also obtain constraints on the couplings cy h R and cs L R from x or y for different 
d u , allowing unparticle contributions to saturate the experimental upper bound on or 
yo- As the contributions from unparticle are suppressed by factors of {m 2 D / A. u ) dhl ~' 1 and 
(m 2 D / A u ) du for vector and scalar unparticle respectively, if n is large the contributions are 
negligible for a fixed Au- We will therefore just consider the lowest possibilities with phase 
7rdu covering all four quadrants plotting results for du in the range of 1 to 3 for illustration 
in Fig. [2] (solid curves) with An fixed to be 1 TeV. In this range, the phase ndu will cover all 
four quadrants. At du equal to half integers (1.5, 2.5), there is no contribution to Am 11 , and 
therefore there are no constraints on cy,s- This is indicated by the two peaks at du = 1-5 and 
2.5 in Fig. [2j At du equal to 1, sin(7r<i w ) = 0, naively the contribution blows off. However, 
at du = 1, A-d u l sin(7TG^) is finite and therefore there is finite contribution as should be for 
a dimension one particle. For other integers, the contribution blows off. For these reasons, 
when reading Fig. [21 one should not taken values too close to integers larger than 1 and half 
integers for du- One can also use data on y^ to constrain cy,s. The results are shown in 
Fig. [3] (solid curve). Note that in this case, at du = 1, there is no contribution to yr> because 
A du = 0. 

From Figs. [2] and [3] (solid curves), we see that constraints from xn and yn give similar 
constraints for cy.s when away from integers and half-integers. The constraint for cs is 
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weaker than cy because of the relative suppression factor vr? D jK u as pointed out earlier. 
Smaller du give stronger constraints on cy,s- 
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FIG. 2: Constraints on the coupling cs (upper curves) and cy (lower curves) as functions of du 
from meson-antimeson oscillation with Ay = 1 TeV by fitting the mass differences Am. The solid 
(blue), dashed (red), and dotted (black) curves are for D° — Z)°, Bd — B& and B s — B S: respectively. 



B° d - B G d System 

Am^ has been measured to be (0.507 ± 0.005)p S - 1 (x = 0.776 ± 0.008) Q, Q- A 
non-zero Ar# d (c/) has not been established. The SM prediction for Am^ agrees with data 



very well. The prediction for is — (26. 7^5) x 10 _4 ps _1 [14| which is very small to be 

measured experimentally. Since SM prediction for Arn^ agrees with data, in our analysis 
to obtain constraints on the couplings cy,s we will just allow the unparticle contributions to 
vary within 2a of experimental error bar. 

Since the SM prediction for Mj^ is positive, if the unparticle contribution is required 
to increase (decrease) the value for Am^ relative to the SM value, then ndu needs to 
be in second and fourth (first and third) quadrants. In Fig. [21 we show the results for 
constraints on cy,s assuming the unparticle contributions saturate 2cr of experimental error 
bar for Am^. The constraints are similar to those obtained from D° — D° oscillation. The 
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FIG. 3: Constraints on the coupling cs (upper curves) and cy (lower curves) from meson-antimeson 
oscillation with An = 1 TeV by fitting the width differences Ar. The solid (blue) and dashed (red) 
curves are for D° — D° and B s — B s , respectively. 

predicted value for y^ d is given by y^ d = x^ d tan(7rt4/) which can be as large as present 
experimental upper bound since at du close to half integers cot(ndu) can be very large, and 
at half integers there is no constraint from Am. Future experiments may tell us more. 
B° s - B° s System 

For B° s - B° s , Ara» is measured to be (17.77 ± 0.12)ps" 1 Q, and AT Bs = 
-(O.OS^^Dps" 1 UJlJ. SM best fits are Am Bs = (19.3 ± 6.68)ps' 1 and Ar Bs = 
— (0.096 ± 0.039)ps _1 [14]. There are differences for central values of SM predictions and 
experimental measurements. Taking these central values and attributing the differences are 
due to unparticle effects, one would favor ndy to be in the first or third quadrants. Since 
both SM predictions and experimental measurements have large errors and they agree with 
in error bars, we will present our constraints on cy,s taking, again, 2a experimental error 
bars for both Am and AT. The results are shown in Figs. [2] and [31 The constraints on cy,s 
are similar to those obtained from D° — D° oscillation. 
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III. MUONIUM-ANTIMUONIUM OSCILLATION 



For muonium and antimuonium oscillation to occur, there must be flavor changing inter- 
actions. To the lowest order, the following unparticle and particle interaction operators will 
contribute, 

a) : \> LL Alr d "L L ^L L 0^ \' EE klr du E R i»E R Oh 

b) : i\ LL A u d »L L ^D»L L O u , i\ EE A u d «E Rl ^E R O u - 

c) : i~\ LL Au du L L ^L L d»O u , i~X EE A^E RllM E R d^O u ; 

d) : \ YE A^L L HE R O u . (11) 

At tree level, exchange of unparticles will generate /zl^e/il^e type of matrix elements. 
The operators has the same form given in eq.(EI) with appropriate replacements of quarks 
by letptons and the associated couplings. 

The SM prediction for muonium and antimuonium oscillation is extremely small. Obser- 
vation of this oscillation at a substantially larger rate will be an indication of new physics. 
Experimentally, no oscillation has been observed. The current upper limit for the probability 
of spontaneous muonium to antimuonium conversion was established at Pg M < 8.3 x 10~ n 
(90% C.L.) in 0.1 T magnetic field Q. 

In the absence of external electromagnetic fields, the probability Pmm °f observing 
a transition can be written as Q Pm M (0T) ~ \5\ 2 /(2Tl), where 5 = 2(M\H eff \M) 
and r M is the muon decay width. Here the effective Hamiltonian is defined as H e ff = 
(Gmm / V2)pT 1 ejiT2e. For the Ti x T 2 = (V ± A) 2 type Hamiltonian (//7a (1 — 7s)e) 2 , the 
transition amplitude is given by 5 = 16G^ M /(v^2va 3 ) for both triplet and singlet muo- 
nium states, where a ~ (am e ) _1 is the Bohr radius. But for (S ± P) 2 type we have 
5 = — 4G A 7 /M /(v / 27ra 3 ) for both triplet and singlet muonium jis| . 

As for our case, omitting m e , the contributions corresponding to parameters Cy hR are 
(V ± A) 2 + (S ± P) 2 type and to parameters cs L R are (S ± P) 2 type. Therefore we have 6 
given by 5 = 12GM M /(^/2na 3 ) and 5 = — 4G A f M /(v^2vra 3 ) for cy L R and cs LR respectively. 

It is important to note that the probability Pmm has strong magnetic field dependence 
which usually occurs in experimental situation. With an external magnetic field, there is a 
reduction factor S R , i-e. P MM (B) = SbPmm(^)- The magnetic field correction factor S R 
describes the suppression of the conversion in the external magnetic field due to the removal 



of degeneracy between corresponding levels in M and M. One has Sr = 0.35 for (V ± A) 2 

n n 

and (S ± P) 2 type interactions at B — 0.1T [16j, |19|]. Using this experimental information, 
one obtains the usual constraint G M m < 3.0 x 10~ 3 Gp for (V ± A) 2 type interaction [16 ]. 
Applying to our case we can put constraints on the relevant parameters and obtain 

A du fmlY"- 1 Me N2| ^ 4.0 x 10- 3 G F 



& r < 



Wsm(d U 7T)m 2 M \A 2 U ) v l,r) — ^ 



A >>u (<iY%^e , 2 , < 1-2 X 10~ 2 G F ^ (r)) 



l 16sm(d U 7r)m 2 M \A 2 U ) [ ° Sl ' rJ 1 " y/2 

where Gp is the Fermi constant. 

Using eq. fU2l) . one can obtain constraints on cy t s f° r given Ky and du- The constraints for 
cy,s are shown in FigJUfor Ay = 1 TeV. At du equal to integers larger than 1 the contribution 
to 5 blows off due to the appearance of sm(irdu) in the denominator of eq. ([!)]) and therefore 
one should take values away from du close to integers. We see that stringent constraint can 
be obtained on cy for small du- The constraint for cs is weak because the suppression factor 
of m^/Aif compared with that for Cy. In general the constraints are weaker compared with 
those obtained from meson and antimeson oscillations since the suppression factors are now 
(m 2 l /A' u ) du ~ 1 and (m 2 l /A u ) du for vector and scalar unparticle contributions which are more 
severe than that for meson-antimeson oscillation cases. 



IV. DISCUSSIONS AND CONCLUSIONS 

In our previous discussions, we have assumed that there is no CP violation in the inter- 
actions between unparticles and particles. We now briefly comment on some implications 
for CP violation. If the parameters A, A and Ay are complex, CP is violated. There may be 
chance to have large enough CP asymmetry for A$l — (1 ~ |<?/p| 4 )/(l + Iq/pI 4 ) which may 
be observed by measuring meson and antimeson semi-leptonic decays. Let us take — B& 
mixing for discussions since there are a large number of B& mesons produced at B-factories 
at KEK and SLAC and more detailed study could be carried out in the near future. 

In the SM, Asl for Bd — Bd system is predicted to be very small (< 10~ 3 ). The reasons 
for this are two folds: small and small relative CP violating phase between M 12 ad 
Ti2. With unparticle interactions, can be much larger than the SM prediction as shown 
before, and with non-zero CP violating phases for A, A and Ay relative CP violating phase 



10 



O 




FIG. 4: Constraints on the coupling cy,s from muonium-antimuonium oscillation with Ay = 1 
TeV. The upper and lower curves are for cs and cy respectively. 

between Ami2 and Ti 2 can be generated. It is possible to have a sizeable Asx- To a good 
approximation, we have, Asl ~ |r^ 2 | sin 0/ 1 M*!f \ . Here is a relative CP violating phase 

between total M\% tal and r^ 9 which is unknown. 

□ 

Experimentally, Asl is constrained to be [12[ —0.0049 ± 0.0038. Allowing the unparticle 
contribution to saturate experimental upper bound on yB d , Asl can easily reach present 
constraint. Similar situation occurs for CP violation in B s — B s system. Measurements of 
Asl can also provide information about unparticle interactions. 

To summarize, we have studied unparticle effects on particle and antiparticle oscillations 
in meson-antimeson, and muonium-antimuonium systems. We found that unlike usual tree 
level contributions to meson oscillations from heavy particle exchange with small T^, the 
unparticle may have sizeable contributions to both Myi and Ti2 due to the fractional di- 
mension du of the unparticle. Numerically we found that very stringent constraints on the 
unparticle and particle interactions can be obtained. If unparticle effect dominates the con- 
tributions (which may happen in D° — D° mixing) to meson mixing parameters x and y, 
x/y = cotijidu)- New constraints on unparticle and particle interactions can also be ob- 



11 



tained using muonium and antimuonium oscillation data. Unparticle interactions can also 
induce large CP violation in meson oscillations. 
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